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Abstract 



Most state-of-the-art techniques for multi-class image segmentation and labeling 
use conditional random fields defined over pixels or image regions. While region- 
level models often feature dense pairwise connectivity, pixel-level models are con- 
siderably larger and have only permitted sparse graph structures. In this paper, we 
consider fully connected CRF models defined on the complete set of pixels in an 
image. The resulting graphs have billions of edges, making traditional inference 
algorithms impractical. Our main contribution is a highly efficient approximate 
inference algorithm for fully connected CRF models in which the pairwise edge 
potentials are defined by a linear combination of Gaussian kernels. Our experi- 
ments demonstrate that dense connectivity at the pixel level substantially improves 
segmentation and labeling accuracy. 



1 Introduction 

Multi-class image segmentation and labeling is one of the most challenging and actively studied 
problems in computer vision. The goal is to label every pixel in the image with one of several prede- 
termined object categories, thus concurrently performing recognition and segmentation of multiple 
object classes. A common approach is to pose this problem as maximum a posteriori (MAP) infer- 
ence in a conditional random field (CRF) defined over pixels or image patches [8, 12, 18, 19, 9]. 
The CRF potentials incorporate smoothness terms that maxinuze label agreement between similar 
pixels, and can integrate more elaborate terms that model contextual relationships between object 
classes. 

Basic CRF models are composed of unary potentials on individual pixels or image patches and pair- 
wise potentials on neighboring pixels or patches [19, 23, 7, 5]. The resulting adjacency CRF struc- 
ture is limited in its ability to model long-range connections within the image and generally results 
in excessive smoothing of object boundaries. In order to improve segmentation and labeling accu- 
racy, researchers have expanded the basic CRF framework to incorporate hierarchical connectivity 
and higher-order potentials defined on image regions [8, 12, 9, 13]. However, the accuracy of these 
approaches is necessarily restricted by the accuracy of unsupervised image segmentation, which is 
used to compute the regions on which the model operates. This limits the ability of region-based 
approaches to produce accurate label assignments around complex object boundaries, although sig- 
nificant progress has been made [9, 13, 14]. 

In this paper, we explore a different model structure for accurate semantic segmentation and labeling. 
We use 8i fully connected CRF that establishes pairwise potentials on all pairs of pixels in the image. 
Fully connected CRFs have been used for semantic image labeling in the past [18, 22, 6, 17], but the 
complexity of inference in fully connected models has restricted their application to sets of hundreds 
of image regions or fewer. The segmentation accuracy achieved by these approaches is again limited 
by the unsupervised segmentation that produces the regions. In contrast, our model connects all 
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MCMC inference, 36 hrs our approach, 0.2 seconds 



Figure 1: Pixel-level classification with a fully connected CRF. (a) Input image from the MSRC-21 dataset. (b) 
The response of unary classifiers used by our models, (c) Classification produced by the Robust CRF [9]. 
(d) Classification produced by MCMC inference [17] in a fully connected pixel-level CRF model; the algorithm 
was run for 36 hours and only partially converged for the bottom image, (e) Classification produced by our 
inference algorithm in the fully connected model in 0.2 seconds. 

pairs of individual pixels in the image, enabling greatly refined segmentation and labeling. The 
main challenge is the size of the model, which has tens of thousands of nodes and billions of edges 
even on low-resolution images. 

Our main contribution is a highly efficient inference algorithm for fully connected CRF models in 
which the pairwise edge potentials are defined by a linear combination of Gaussian kernels in an ar- 
bitrary feature space. The algorithm is based on a mean field approximation to the CRF distribution. 
This approximation is iteratively optimized through a series of message passing steps, each of which 
updates a single variable by aggregating information from all other variables. We show that a mean 
field update of all variables in a fully connected CRF can be performed using Gaussian filtering 
in feature space. This allows us to reduce the computational complexity of message passing from 
quadratic to linear in the number of variables by employing efficient approximate high-dimensional 
filtering [16, 2, 1]. The resulting approximate inference algorithm is sublinear in the number of 
edges in the model. 

Figure 1 demonstrates the benefits of the presented algorithm on two images from the MSRC-21 
dataset for multi-class image segmentation and labeling. Figure 1(d) shows the results of approxi- 
mate MCMC inference in fully connected CRFs on these images [17]. The MCMC procedure was 
run for 36 hours and only partially converged for the bottom image. We have also experimented with 
graph cut inference in the fully connected models [11], but it did not converge within 72 hours. In 
contrast, a single-threaded implementation of our algorithm produces a detailed pixel-level labeling 
in 0.2 seconds, as shown in Figure 1(e). A quantitative evaluation on the MSRC-21 and the PAS- 
CAL VOC 2010 datasets is provided in Section 6. To the best of our knowledge, we are the first to 
demonstrate efficient inference in fully connected CRF models at the pixel level. 

2 The Fully Connected CRF Model 

Consider a random field X defined over a set of variables {Xi, . . . , X^}- The domain of each 
variable is a set of labels >C = {/i, ^2, • • • , ^/c}- Consider also a random field I defined over variables 
{/i, . . . , /at}. In our setting, I ranges over possible input images of size N and X ranges over 
possible pixel-level image labelings. Ij is the color vector of pixel j and Xj is the label assigned to 
pixel j. 

A conditional random field (I, X) is characterized by a Gibbs distribution 
P(X|I) = wiyexp(— ^^^^ 0c(Xc|I)), where Q = (V, f) is a graph on X and each clique c 
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in a set of cliques Cg in Q induces a potential (pc [15]. The Gibbs energy of a labeling x G 
is ^(x|I) = Y2ceCg <^c(xc|I). The maximum a posteriori (MAP) labeling of the random field is 
X* = arg maxxe/:^^ P(x|I). For notational convenience we will omit the conditioning in the rest of 
the paper and use ipd'^c) to denote (/)c(xc|I). 

In the fully connected pairwise CRF model, Q is the complete graph on X and Cg is the set of all 
unary and pairwise cliques. The corresponding Gibbs energy is 

^(x) = ^^iljp{xi,Xj), (1) 

i i<j 

where i and j range from 1 to N. The unary potential ipu{xi) is computed independently for each 
pixel by a classifier that produces a distribution over the label assignment Xi given image features. 
The unary potential used in our implementation incorporates shape, texture, location, and color 
descriptors and is described in Section 5. Since the output of the unary classifier for each pixel 
is produced independently from the outputs of the classifiers for other pixels, the MAP labeling 
produced by the unary classifiers alone is generally noisy and inconsistent, as shown in Figure 1(b). 

The pairwise potentials in our model have the form 

i;j,{xi,Xj) = f,{xi,Xj)j:^^,w^^^k(^\U,fj), (2) 

V ' 

/C(f,,f,) 

where each /c*^^^ is a Gaussian kernel /c^^^ (f^, fj) = exp(— ^ (f^ — fjj^A*^"^^ (f^ — fj)), the vectors 
and fj are feature vectors for pixels i and j in an arbitrary feature space, w^^"^ are linear combination 
weights, and /i is a label compatibility function. Each kernel k^^"^ is characterized by a symmetric, 
positive-definite precision matrix A^'^^ which defines its shape. 

For multi-class image segmentation and labeling we use contrast- sensitive two-kernel potentials, 
defined in terms of the color vectors li and Ij and positions pi and pj : 

appearance kernel smoothness kernel 

The appearance kernel is inspired by the observation that nearby pixels with similar color are likely 
to be in the same class. The degrees of nearness and similarity are controlled by parameters Qq, and 
Q^. The smoothness kernel removes small isolated regions [19]. The parameters are learned from 
data, as described in Section 4. 

A simple label compatibility function /i is given by the Potts model, fi{xi,Xj) = [xi ^ Xj]. It 
introduces a penalty for nearby similar pixels that are assigned different labels. While this simple 
model works well in practice, it is insensitive to compatibility between labels. For example, it 
penalizes a pair of nearby pixels labeled "sky" and "bird" to the same extent as pixels labeled "sky" 
and "cat". We can instead learn a general symmetric compatibility function fi{xi^Xj) that takes 
interactions between labels into account, as described in Section 4. 

3 Efficient Inference in Fully Connected CRFs 

Our algorithm is based on a mean field approximation to the CRF distribution. This approxima- 
tion yields an iterative message passing algorithm for approximate inference. Our key observation 
is that message passing in the presented model can be performed using Gaussian filtering in fea- 
ture space. This enables us to utilize highly efficient approximations for high-dimensional filtering, 
which reduce the complexity of message passing from quadratic to linear, resulting in an approxi- 
mate inference algorithm for fully connected CRFs that is linear in the number of variables and 
sublinear in the number of edges in the model. 

3.1 Mean Field Approximation 

Instead of computing the exact distribution P(X), the mean field approximation computes a dis- 
tribution Q(X) that minimizes the KL-divergence 'D{Q\\P) among all distributions Q that can be 
expressed as a product of independent marginals, Q(X) = Yli Qi{Xi) [10]. 
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Minimizing the KL-divergence, while constraining Q(X) and Qi{Xi) to be vahd distributions, 
yields the following iterative update equation: 

Qii^i = = ^exp I -M^,) - ^ ^^W(f,,f,)o,(r) i . (4) 

* [ I'eC m=l j^i J 

A detailed derivation of Equation 4 is given in the supplementary material. This update equation 
leads to the following inference algorithm: 

Algorithm 1 Mean field in fully connected CRFs 

> Qi{xi) ^ ^ exp{-(t)u{xi)} 

> See Section 6 for convergence analysis 

> Message passing from all Xj to all Xi 
> Compatibility transform 

> Local update 



Each iteration of Algorithm 1 performs a message passing step, a compatibility transform, and a 
local update. Both the compatibility transform and the local update run in linear time and are highly 
efficient. The computational bottleneck is message passing. For each variable, this step requires 
evaluating a sum over all other variables. A naive implementation thus has quadratic complexity in 
the number of variables N. Next, we show how approximate high-dimensional filtering can be used 
to reduce the computational cost of message passing to linear. 

3.2 Efficient Message Passing Using High-Dimensional Filtering 

From a signal processing standpoint, the message passing step can be expressed as a convolution 
with a Gaussian kernel Gj^(m,) in feature space: 

Qi"\i) = E,-ev^^"^Hf.,f.)Q.(0 - Qiii) = [Ga(^) ^ Q(0] (fi) -Qiii)- (5) 

message passing Qi^^ (I) 

We subtract Qi{l) from the convolved function Q-"^^ (/) because the convolution sums over all vari- 
ables, while message passing does not sum over Qi. 

This convolution performs a low-pass filter, essentially band-limiting q[^\i). By the sampling 
theorem, this function can be reconstructed from a set of samples whose spacing is proportional 
to the standard deviation of the filter [20]. We can thus perform the convolution by downsampling 
Q{1), convolving the samples with Ga(^) » upsampling the result at the feature points [16]. 



Algorithm 2 Efficient message passing: Q (/) = X^^^y A^*^^^ 


'(fi,fi)Qi(0 


Ql{l) ^ downsample((5(/)) 
Q^"'\l) ^ upsample(oj."^(0) 


> Downsample 
> Convolution on samples f^^ 
> Upsample 



A conmion approximation to the Gaussian kernel is a truncated Gaussian, where all values beyond 
two standard deviations are set to zero. Since the spacing of the samples is proportional to the stan- 
dard deviation, the support of the truncated kernel contains only a constant number of sample points. 
Thus the convolution can be approximately computed at each sample by aggregating values from 
only a constant number of neighboring samples. This implies that approximate message passing can 
be performed in 0{N) time [16]. 

High-dimensional filtering algorithms that follow this approach can still have computational com- 
plexity exponential in d. However, a clever filtering scheme can reduce the complexity of the con- 
volution operation to 0{Nd). We use the permutohedral lattice, a highly efficient convolution data 



Initialize Q 

while not converged do 

Qt\l) ^ E,-/i^^"^nfi,fi)Oi(0 for all m 

Qi{xi) ^ exp{-V^^(x^) - Qi{xi)} 
normalize Qi{xi) 
end while 
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structure that tiles the feature space with simpHces arranged along d-\-l axes [1]. The permutohedral 
lattice exploits the separability of unit variance Gaussian kernels. Thus we need to apply a whitening 
transform f = [/f to the feature space in order to use it. The whitening transformation is found us- 
ing the Cholesky decomposition of A^^^ into UU^ . In the transformed space, the high-dimensional 
convolution can be separated into a sequence of one-dimensional convolutions along the axes of the 
lattice. The resulting approximate message passing procedure is highly efficient even with a fully 
sequential implementation that does not make use of parallelism or the streaming capabilities of 
graphics hardware, which can provide further acceleration if desired. 



4 Learning 

We learn the parameters of the model by piecewise training. First, the boosted unary classifiers are 
trained using the JointBoost algorithm [21], using the features described in Section 5. Next we learn 
the appearance kernel parameters w^^\ 0^, and Ojs for the Potts model, w^^^ can be found efficiently 
by a combination of expectation maximization and high-dimensional filtering. Unfortunately, the 
kernel widths 6 a and cannot be computed effectively with this approach, since their gradient 
involves a sum of non-Gaussian kernels, which are not amenable to the same acceleration techniques. 
We found it to be more efficient to use grid search on a holdout validation set for all three kernel 
parameters w^^\Oa and ^^3. 

The smoothness kernel parameters w^'^'^ and Oj do not significantly affect classification accuracy, 
but yield a small visual improvement. We found w^'^'^ = 0^ = 1 to work well in practice. 

The compatibility parameters /i(a, b) = a) are learned using L-BFGS to maximize the log- 
likelihood ^(/i : X, T) of the model for a validation set of images X with corresponding ground 
truth labelings T. L-BFGS requires the computation of the gradient of £, which is intractable to 
estimate exactly, since it requires computing the gradient of the partition function Z. Instead, we 
use the mean field approximation described in Section 3 to estimate the gradient of Z. This leads to 
a simple approximation of the gradient for each training image: 

(6) 

where (X(^),r(^)) is a single training image with its ground truth labeling and T*^^^ (a) is a binary 

image in which the ith pixel T^^'' (a) has value 1 if the ground truth label at the ith pixel of T^^^ is 
a and otherwise. A detailed derivation of Equation 6 is given in the supplementary material. 

The sums ^j^i ij)Tj{b) and ^i^i^ ^j)Qi{^) t)oth computationally expensive to eval- 
uate directly. As in Section 3.2, we use high-dimensional filtering to compute both sums efficiently. 
The runtime of the final learning algorithm is linear in the number of variables TV. 



5 Implementation 

The unary potentials used in our implementation are derived from TextonBoost [19, 13]. We use 
the 17 -dimensional filter bank suggested by Shotton et al. [19], and follow Ladicky et al. [13] by 
adding color, histogram of oriented gradients (HOG), and pixel location features. Our evaluation 
on the MSRC-21 dataset uses this extended version of TextonBoost for the unary potentials. For 
the VOC 2010 dataset we include the response of bounding box object detectors [4] for each object 
class as 20 additional features. This increases the performance of the unary classifiers on the VOC 
2010 from 13% to 22%. We gain an additional 5% by training a logistic regression classifier on the 
responses of the boosted classifier. 

For efficient high-dimensional filtering, we use a publicly available implementation of the permuto- 
hedral lattice [1]. We found a downsampling rate of one standard deviation to work best for all our 
experiments. Sampling-based filtering algorithms underestimate the edge strength A:(fj, fj) for very 
similar feature points. Proper normalization can cancel out most of this error. The permutohedral 
lattice allows for two types of normalizations. A global normalization by the average kernel strength 
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k = j k{fi^ fj) can correct for constant error. A pixelwise normalization by ki = k{fi^ fj) 

handles regional errors as well, but violates the CRF symmetry assumption ipp {xi , ) = ipp {xj , x^) . 
We found the pixelwise normalization to work better in practice. 



6 Evaluation 

We evaluate the presented algorithm on two standard benchmarks for multi-class image segmen- 
tation and labeling. The first is the MSRC-21 dataset, which consists of 591 color images of size 
320 X 213 with corresponding ground truth labelings of 21 object classes [19]. The second is the 
PASCAL VOC 2010 dataset, which contains 1928 color images of size approximately 500 x 400, 
with a total of 20 object classes and one background class [3]. The presented approach was evalu- 
ated alongside the adjacency (grid) CRF of Shotton et al. [19] and the Robust CRF of KohH et 
al. [9], using publicly available reference implementations. To ensure a fair comparison, all models 
used the unary potentials described in Section 5. All experiments were conducted on an Intel 17-930 
processor clocked at 2.80GHz. Eight CPU cores were used for training; all other experiments were 
performed on a single core. The inference algorithm was implemented in a single CPU thread. 

Convergence. We first evaluate the convergence of the mean field approximation by analyzing 
the KL-divergence between Q and P. Figure 2 shows the KL-divergence between Q and P over 
successive iterations of the inference algorithm. The KL-divergence was estimated up to a constant 
as described in supplementary material. Results are shown for different standard deviations Oa and 
0f3 of the kernels. The graphs were aligned at 20 iterations for visual comparison. The number of 
iterations was set to 10 in all subsequent experiments. 

MSRC-21 dataset. We use the standard split of the dataset into 45% training, 10% vaHdation and 
45% test images [19]. The unary potentials were learned on the training set, while the parameters of 
all CRF models were learned using holdout validation. The total CRF training time was 40 minutes. 
The learned label compatibility function performed on par with the Potts model on this dataset. 
Figure 3 provides qualitative and quantitative results on the dataset. We report the standard measures 
of multi-class segmentation accuracy: "global" denotes the overall percentage of correctly classified 
image pixels and "average" is the unweighted average of per-category classification accuracy [19, 9]. 
The presented inference algorithm on the fully connected CRF significantly outperforms the other 
models, evaluated against the standard ground truth data provided with the dataset. 

The ground truth labelings provided with the MSRC-21 dataset are quite imprecise. In particular, 
regions around object boundaries are often left unlabeled. This makes it difficult to quantitatively 
evaluate the performance of algorithms that strive for pixel-level accuracy. Following Kohli et al. [9], 
we manually produced accurate segmentations and labelings for a set of images from the MSRC-21 
dataset. Each image was fully annotated at the pixel level, with careful labeling around complex 
boundaries. This labeling was performed by hand for 94 representative images from the MSRC- 
21 dataset. Labeling a single image took 30 minutes on average. A number of images from this 
"accurate ground truth" set are shown in Figure 3. Figure 3 reports segmentation accuracy against 
this ground truth data alongside the evaluation against the standard ground truth. The results were 
obtained using 5-fold cross validation, where | of the 94 images were used to train the CRF pa- 
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Figure 2: Convergence analysis, (a) KL-divergence of the mean field approximation during successive itera- 
tions of the inference algorithm, averaged across 94 images from the MSRC-21 dataset. (b) Visualization of 
convergence on distributions for two class labels over an image from the dataset. 
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Figure 3: Qualitative and quantitative results on the MSRC-21 dataset. 

rameters. The unary potentials were learned on a separate training set that did not include the 94 
accurately annotated images. 

We also adopt the methodology proposed by Kohli et al. [9] for evaluating segmentation accuracy 
around boundaries. Specifically, we count the relative number of misclassified pixels within a nar- 
row band ("trimap") surrounding actual object boundaries, obtained from the accurate ground truth 
images. As shown in Figure 4, our algorithm outperforms previous work across all trimap widths. 

PASCAL VOC 2010. Due to the lack of a publicly available ground truth labeling for the test 
set in the PASCAL VOC 2010, we use the training and validation data for all our experiments. We 
randomly partitioned the images into 3 groups: 40% training, 15% validation, and 45% test set. Seg- 
mentation accuracy was measured using the standard VOC measure [3]. The unary potentials were 
learned on the training set and yielded an average classification accuracy of 27.6%. The parameters 
for the Potts potentials in the fully connected CRF model were learned on the validation set. The 
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Figure 4: Segmentation accuracy around object boundaries, (a) Visualization of the "trimap" measure, (b) 
Percent of misclassified pixels within trimaps of different widths. 
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Figure 5: Qualitative results on the PASCAL VOC 2010 dataset. Average segmentation accuracy was 30.2%. 

fully connected model with Potts potentials yielded an average classification accuracy of 29.1%. 
The label compatibility function, learned on the validation set, further increased the classification 
accuracy to 30.2%. For comparison, the grid CRF achieves 28.3%. Training time was 2.5 hours and 
inference time is 0.5 seconds. Qualitative results are provided in Figure 5. 

Long-range connections. We have examined the value of long-range connections in our model by 
varying the spatial and color ranges 6 a and 6f3 of the appearance kernel and analyzing the resulting 
classification accuracy. For this experiment, w'^^^ was held constant and w'^'^^ was set to 0. The 
results are shown in Figure 6. Accuracy steadily increases as longer-range connections are added, 
peaking at spatial standard deviation ofOa = Ql pixels and color standard deviation 6>/3 = 11. At this 
setting, more than 50% of the pairwise potential energy in the model was assigned to edges of length 
35 pixels or higher. However, long-range connections can also propagate misleading information, 
as shown in Figure 7. 




(a) Quantitative (b) Qualitative 

Figure 6: Influence of long-range connections on classification accuracy, (a) Global classification accuracy on 
the 94 MSRC images with accurate ground truth, as a function of kernel parameters Oa and O13. (b) Results for 
one image across two slices in parameter space, shown as black lines in (a). 



Discussion. We have presented a highly efficient approximate inference algorithm for fully con- 
nected CRF models. Our results demonstrate that dense pixel-level connectivity leads to signif- 
icantly more accurate pixel-level classification performance. Our single-threaded implementation 
processes benchmark images in a fraction of a second and the algorithm can be parallelized for 
further performance gains. 
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Figure 7: Failure cases on images from the PASCAL VOC 2010 (left) and the MSRC-21 (right). Long-range 
connections propagated misleading information, eroding the bird wing in the left image and corrupting the legs 
of the cat on the right. 
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